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Abstract 



Results of a large scale numerical simulation show that the high temperature Chern-Simons 
number diffusion rate in the electroweak theory has a classical limit T = ^(a^r)^, where 
K = 1.09 lb 0.04 and a,„ is the weak fine structure constant. 
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1 Introduction 



Topology-changing transitions play an important role in the electroweak theory since they 
are accompanied by baryon-number violating processes [Q. At zero temperature these 
processes can only appear as true quantum effects since the energy of the sphaleron con- 
figuration is a barrier separating two gauge equivalent vacua ||2|. At temperatures far above 
the sphaleron energy one would expect that this barrier height, as well as the concept of 
vacuum to vacuum amplitude, are irrelevant for the transition rate. Unfortunately there 
presently is no way to calculate analytically the transition rate at these high temperatures. 
From a numerical point of view the situation is not any better if we insist on a fully quan- 
tum treatment, since the only non-perturbative tool available is lattice gauge theories in 
the imaginary time formalism. Real time transition rates are not easily extracted in this 
formalism. 

The object of our study is the topological charge 

= 3^i7 ^'^^;^^"'^'^ = ^cs(t)-A^cs(0), (1) 
where the Chern-Simons number 

Ncsit) ^ ^ e,,k J d'x (^F^^At - ^e'^'^A^A'^Al^ . (2) 

The theory has a periodic (with period 1) potential in the A'cs direction. This periodicity 
is due to the symmetry under large gauge transformations, changing A'cs by an integer. 
There also are dynamical processes whose result for a field configuration is the same as a 
large gauge transformation. These processes have an integer topological charge. Due to 
nonlinear interactions of A'cs with other degrees of freedom of the theory it is expected 
that for large times t B{t) is a random walk in a periodic A'cs potential 

{B\t))^ = rvt, (3) 

where V is the space volume, by ()y we mean average over the canonical ensemble, and 
r, the diffusion constant per unit volume, is the transition rate we are interested in Q. 

At high temperatures processes with an integer topological charge occur predominantly 
by classically allowed thermal activation rather than by quantum tunneling. Some time 
ago it was suggested to exploit this property and to determine F entirely within the 
classical theory. This is done by solving classical equations of motion for the fields and 
averaging the resulting B'^{t) over the classical canonical ensemble. Behind this approach 
is the expectation that the topology-changing transitions mostly involve field fluctuations 
with long wavelengths and large magnitudes (hence many quanta). Such fluctuations are 
well described by classical statistical mechanics. An argument usually given to support 
this picture is as follows. Consider a classically allowed process with B{t) = 1, whereby 
the field configuration evolves from the vicinity of one minimum in the N^s potential to 
the vicinity of another one. If a typical linear space-time size of the process is r, the 
system will cross an energy barrier of height E ~ l/ar, where a = /^tt. While large r 
are favored energetically, at high temperatures r is unlikely to exceed the inverse of the 
magnetic screening mass mmag ~ olT. We then expect that processes with r ~ 1/aT give 
the most important contribution to F. In the electroweak theory, where « 1/30, these 
r are much larger than the thermal wavelength, i.e., fluctuations of size r can be treated 
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classically. In order to estimate the transition rate, the Boltzmann factor associated with 
the barrier, exp(— i?/T) ~ 1, is divided by the space-time volume of the process, to 
give [| 

r = K{aT)^, (4) 

where k is a dimensionless constant. This estimate should be valid at temperatures far 
above the electroweak transition scale. We also expect that at such high temperatures the 
mass scales of the Higgs sector are irrelevant, and that T can be reliably obtained from a 
pure Yang-Mills theory. 

The classical approximation offers an enormous simplification in computing real-time 
correlation functions. However, it is well known that classical field theory at finite temper- 
ature suffers from the Rayleigh- Jeans ultraviolet divergence. This means that, upon cutoff 
regularization, classical thermodynamic average of a generic observable will be sensitive 
to thermal fluctuations at the cutoff scale A Q. The classical rate F, if it is to be trusted, 
must not show such sensitivity. This is indeed what we expect: considerations leading 
to the estimate (^) may be repeated for a purely classical theory which also possesses a 
magnetic mass mmag ~ oT. In the quantum case, we required separation between the 
magnetic and the thermal scales: qT <^ T. Likewise, in the classical theory we should 
require separation between the magnetic and the cutoff scales: aT <^ A. These two con- 
ditions make it plausible that T is not strongly affected by either quantum corrections or 
the classical cutoff. 

On the other hand, it is not clear how the T independence of A is technically possible: 
after all, B{t), the topological charge per unit time, is an ultraviolet divergent quantity 
whose classical standard deviation is given perturbatively by 

iB?)^=(^y (( [d'xEnx)BnA{x))y\ -rW, (5) 
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where Ef and Bf are, respectively, color electric and magnetic fields, and V is the space 
volume. One then naturally expects the Rayleigh-Jeans divergency to show up in (^B'^{t))rp 
itself. 

The key point, confirmed by the numerical simulations in the following, is that the UV 
divergency shows up in {B'^{t))j, in a way that does not affect F. To see how this could be 
anticipated, consider the structure of {B'^{t))rp in some more detail. The classical motion 
of A^cs can be thought of as consisting of two pieces: "true" thermal fluctuations and a 
random walk between gauge equivalent vacua. If the theory did not have an inflnite series 
of degenerate vacua, Ncs would simply fluctuate around zero, as it does in an Abelian 
theory. The properties of these thermal fluctuations can be studied perturbatively. As 
discussed in the Appendix, in an Abelian theory, where no random walk of A'^cs is possible, 
{B'^{t))j, is bounded 

/ o \ Abclian / r,\ Abelian 

(BHt))^ <2((iVcs)^)^ (6) 
and at large t fluctuates around the average value 

/ o\ Abelian ( CtT)'^ 

y ^ ((A^cs)')^ « 0.00228^^ y, (7) 

diverging, as expected, linearly with the lattice cutoff A = 1/a. In the nonabelian theory 
the random walk of A'cs is superimposed with this thermal background. A reasonable 
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approximation to the nonabelian {B'^{t))j, is then 



B'{t))^ = a. + TVt, (8) 

where c ~ 3, assuming that thermal fluctuations of A'^cs are exactly those for three inde- 
pendent copies of an Abelian theory. 

If we let A — > oo for fixed t the r.h.s. of (|8|) diverges. On the other hand, for a fixed 
cutoff (§) is well approximated by (|3|) for large enough t and allows unambiguously to 
extract F. We will show that numerical simulations, using a classical thermal ensemble 
with a lattice cutoff A = 1/a, point to a F independent of a for sufficiently small a. On 
dimensional grounds this means that (Q) holds classically in the continuum limit, with k 
of @ being a non-perturbative constant of the classical theory. For the gauge group SU (2) 
we find k = 1.09 it 0.04. This is the central result of our work. 

The ultraviolet insensitivity of F extracted from (^) might be linked to the topological 
aspect of the special field tensor combination F^j^F""^^ in non-abelian field theories. Re- 
placing F^j^F"'^'^ by a generic bilinear of field tensors with zero average would yield a UV 
divergent transition rate. We illustrate this point by measuring in addition the quantity 



r]{t) whose diffusion rate is related to the shear viscosity of the gluon plasma [14|: 

7]{t)= J^dtJd''x{Tu-T22), (9) 
where Tij are components of the energy-momentum tensor. 



2 The model and the numerical method 

The philosophy behind the classical measurements of real-time quantities has been pre- 
sented in detail in a number of articles [^, |8| . Here we only concentrate on the essential 
technical points. We work in the temporal gauge. The classical dynamics of SU(2) Yang- 
Mills theory on the lattice is given by Kogut-Susskind Hamiltonian 

H = \J2 EfEf + E (l - \^Uu) . (10) 

The variable set consists of SU(2) matrices Ui and lattice analogs of electric field Ei, both 
residing on the links of a cubic lattice. We shall also use notation / = j,n for a link 
emanating from site j in positive direction n. The second term in (^) is the standard 
plaquette term, representing the color magnetic energy. Every variable v evolves in time 
according to its canonical equation v = {H,v}, derived using Poisson brackets between 
the independent variables: 

{Ei^,E^,} = 26u'e''^^Ej; {Ei' ,U,} = -i5u'Uia^; {Ui,U,} = (11) 

(no summation over repeated /; a" are the Pauli matrices). The time evolution is subject 
to the set of three Gauss' constraints per site 

E {Eln - Ef^n,r) = 0, (12) 
n 

where 
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These constraints are to be imposed on initial conditions. We use the leapfrog algorithm 
to integrate the equations of motion, performing exact integration of the E equations for 
fixed U and of the U equations for fixed E. In this way the Gauss' laws are obeyed exactly, 
independently of the time step, and the only source of spurious static charge is the limited 
computer accuracy. 

The increment of Chern-Simons number is computed as in ((l|) with the topological 
charge per unit time approximated by 

^cs = ^ E + Ef.n,n) E Tr (^n,,^^") , (13) 

where Oj^^ denotes plaquettes in the plane perpendicular to direction n, running counter- 
clockwise in that plane, and beginning and ending at site j. 

Our purpose is to study the average of B'^{t) {cf (||)) over the thermal ensemble of 
initial field configurations. Due to the presence of first-class constraints (Gauss' laws) 
this point is nontrivial. The initial configurations must satisfy the constraints and be 
thermally distributed in the reduced phase space of physical, gauge-invariant degrees of 
freedom. This goal is achieved by using the constraint-respecting Langevin algorithm. 
The set of Langevin equations reads (/? is the inverse temperature in lattice units) 

Uj,k = {H, Uj,k} + {rj,mn{t) - P{Tj,mn, H}) {Tj^mn, Uj,k}; Ef = {H, E^}, (14) 

where, for fixed j, Tj^mn = \^E'j,m^'j,n (assuming for definitness m < n), and Tj^mnit) is 
a random variable with zero average and correlation 

j ,mnit)'^ j ' ,m' n' {t )) ~ 26(t t )5jji5^Yi,m'n'- 

We are free to choose the friction coefficient 7 > 0. This freedom can be used to optimize 
the algorithm performance. Evolution generated by ( p!4[ ) is not to be confused with real- 
time evolution according to canonical equations of motion: the sole purpose of ( p^ is 
to bring the system in question to thermal equilibrium at temperature The set 

of equations ( [l^ ) is integrated numerically using a modified leapfrog algorithm, which 
obeys the Gauss' laws as accurately as a computer arithmetic allows. For a more detailed 
description of the algorithm we refer to 

Results of simulation must eventually be converted from lattice units to the physical 
ones in the continuum. To do so it is convenient to parametrize the link matrices as 
Ui = exp(iac7"A"/2) where a is the lattice spacing and Af is the gauge potential. It 
follows from the equation of motion for Ui that = iTr(c7"C/;^C//)/2. Therefore for small 
a (|lO|) takes form 

^ « y E {^InKn + "'^"n^^n) , (15) 

where is a chromo-magnetic field. Comparing (|T5| ) to the standard continuum Hamil- 
tonian 

Hc = ^j d\ [AfAf + BtBf) 

we find the correspondence between the time and the temperature expressed in lattice (L) 
and continuum (C) units: 

4 

tc = ati; Tc = -i^Tl. (16) 

Taking into account that A'^cs itself is a dimensionless quantity, (^) on the lattice takes 
form Tl = k{Tl/tt)'^, where the lattice rate is per unit cell. 
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3 Results 



We performed series of simulations using methods described in the previous section for 
three values of the inverse lattice temperature [i = 10, 12, 14 on lattices of equal size 
Ll and periodic boundary conditions in all spatial directions. The value of Ll ranged 
between 8 and 32. Between 20 and 60 statistically independent initial configurations were 
produced for every f3,LL pair (60 for the largest volumes, used for the estimate of k). 
At low temperatures considered the system is only weakly nonlinear, and we expect an 
average energy per physical degree of freedom to be close to 1/2/3. This is indeed what we 
observe: the average chromo-electric energy per site was found to be 3//3 with a three-digit 
accuracy. Our measurement results for the topological charge density are also close to the 
perturbative estimate (p3|). 

Every independent initial configuration was let evolve in real time for 5000 time units. 
We used the integration time step of 0.05 for the classical equations of motion. From our 
experience in earlier work 0] we know that reducing the time step beyond this value has 
negligible impact on the diffusion of Chern-Simons number. During the real-time evolution 
the value of B{t), determined by numerically integrating the lattice topological charge per 
unit time, was recorded with intervals between 0.5 and 1 time unit. In the same way, we 
recorded the time history of 7/(t). Thus our results for {B'^{t))T, discussed in the following, 
represent averaging both over the canonical ensemble of initial configurations and within 
every individual real-time trajectory. The same averaging method was used in earlier work 
on a low-dimensional model . 

A typical plot of {B^it))T 

function of the lag t is shown in Figures |l|and|. As 
expected, for very short lags {t < 1) {B'^{t))T is mostly accounted for by the perturbative 
result obtained from (p^). For larger lags a random walk of A'^cs sets in, and {B'^{t))T 
rapidly approaches a linear asymptote. Thus, for large enough t, {B'^{t))T can be fitted to 
the form (^). The range of t suitable for this linear fit is bounded from below by the initial 
transient behavior, and from above by the growing statistical errors in our determination 
of {B'^{t))T- With the temperatures and lattice sizes used, we found the optimal range to 
be approximately 40 < t < 200. As Figures |l] and |2| indicate, at these lags {B'^{t))T is 
dominated by the nonperturbative effects. The central values of {B'^{t))T for larger t are 
consistent with linear dependence on t. 

We determined the coefficient k of (^) by performing a correlated fit of {B'^{t))T to the 
form (|8|) for data points within the optimal range. As Figure |^ shows, the values of k for 
different combinations of /? and Ll depend, to a good approximation, on a dimensionless 
ratio Ll/(3 (note, in particular, the equality of k values for (3 = 12, Li = 12 and for 
(3 = 14, Ll = 14). We conclude that in the range of temperatures considered k is nearly 
insensitive to the lattice spacing, i.e., close to a finite continuum limit. On the other hand, 
K exhibits strong finite-size effects and grows rapidly for L^/P < 2. This is in agreement 
with our expectation that the diffusion of Acs is dominated by large-scale structures whose 
linear size is of order /3. Once the lattice is too small to accommodate objects of this size, 
the rate rapidly decrease^. For large L^/P > 2 k reaches a constant value. In fact, 
the length scale at which k saturates can be naturally interpreted in terms of magnetic 
mass, nim = 0.466g^T in physical units |11]: k stops growing as soon as the lattice size 

^That the sphaleron-like transitions are related to large scale objects has been verified directly by 
"freezing" the time evolution in the neighborhood of N^s = n+1/2. One finds energy profiles corresponding 
to extended object |loj, in contrast to the situation when one is far from Ncb = n + 1/2. In addition one 
could check that level crossing occurred as expected from continuum physics. 
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Figure 1: Diffusion of Ncs'- mean square of B(t) as a function of the lag t for /? = 12, 
Ll = 32. 
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Figure 2: An enhanced region of < i < 25 of Figure 1. The perturbative estimate of 
{B'^{t))j, is shown by the sohd curve. 
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Figure 3: The rate prefactor k dependence on the dimensionless parameter Li/fi. The 
data points are for /3 = 10 (diamonds), 12 (squares), and 14 (triangles). 



exceeds twice the inverse magnetic mass (the factor of 2 relating the two lengths is to be 
expected with periodic boundary conditions). Our estimate, based on (3 = 12, Li = 32 
is K = 1.09 lb 0.04. Estimates from other large-volume measurements are consistent with 
this one within the error bars. 

Our results for k are in sharp contrast with those for the shear viscosity, extracted 
from the diffusion rate of r/(t). In the chosen range of temperatures r/(i) behaves as 
an underdamped random walk, and determination of the corresponding diffusion constant 
requires knowledge of {ff'{t))T at lags t of the order of 1000 (as compared tot < 200 for the 
Chern-Simons number), making it a difficult task. As Figure ^ shows, the corresponding 
diffusion rate F^ shows no marked temperature or finite-size dependence. The absence 
of finite size effects for F^ indicates that, unlike k, this quantity is not predominantly 
sensitive to long-wavelength modes of the system. Independence of F^ of the temperature 
means, on dimensional grounds, that F^, and hence the classical shear viscosity, diverges 
as a~^, where a is the lattice spacing. This cutoff dependence of transport coefficients is 
found in a similar situation in condensed-matter physics, for phonons at temperatures far 



above the Debye temperature [15 



4 Discussion 

We have seen that k ~ 1 in a classical SU{2) gauge theory, and that we consequently have 
an observable (the diffusion rate between gauge equivalent non-Abelian vacua) which is 
not plagued by the generic ultraviolet thermal fluctuations of classical field theory. In 
addition we expect that the result obtained describes the high temperature limit of the 
classical field theory corresponding to the standard model. 
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Figure 4: Diffusion constant per unit volume of plotted against the dimensionless 
parameter Ll/(3. The data points are for /? = 10 (diamonds), 12 (squares), and 14 
(triangles). 



While this result is interesting by itself, it is of even more interest to understand to 
what extent it is a good approximation to the full quantum theory. In the Introduction we 
gave qualitative arguments in favor of reasonable agreement between the classical result 
for the rate and the full quantum one, essentially reproducing the considerations which 
originally led Grigoriev, Rubakov and Shaposhnikov to suggest the classical real time 
method. Unfortunately, not much is known at present about the size and functional form 
of quantum corrections to dynamical, real-time quantities. The best we can offer is an 
educated guess based on a much better understood relation between the static properties 



of the quantum theory and of the effective dimensinally reduced one fl^, 13 1. The full 
thermal path integral, where the bosonic fields obey periodic boundary conditions in the 
imaginary time direction, reduces in the high temperature limit to a partition function of 
a three-dimensional statistical mechanical problem. This partition function inherits the 
UV thermal cutoff from the underlying quantum field theory. If a lattice regularization is 
used in the quantum theory, the spatial part of this lattice will provide the regularization 
of the thermal three-dimensional theory. It is known that dimensional reduction is not 
completely equivalent to the classical approximation. The correspondence between the 
two is conveniently described in terms of the field component Aq{x) which is the Lagrange 
multiplier of the Gauss' law. The latter may be imposed with the help of a 5-function factor 
in the classical statistical weight. Alternatively, one can introduce Aq{x) and integrate 
out the color electric fields. The classical partition function then takes form 

Z = /pAiP^o e"'^<=«(^''^»\ (17) 
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where 

1 



SesiA) = ^jd'x \F^Ff^ + {DfAlYl . (18) 

On the other hand, the dimensionally reduced partition function has the same form as 
(ITtI), but this time the effective action includes Higgs potential terms for Aq and reads 

SesiA) = I d'x [f-F- + {DfAlf + ^^{T){A-f + cg{Tf{{Alff + • 

m 

where the parameters g{T) and /i(T) are calculable functions of the temperature, g{T) 
being the running, temperature dependent coupling constant, unlike in the purely classical 
case. The important point is that, despite the apparent inequivalence of the two partition 
functions, we do not expect any qualitative difference between the dimensionally reduced 
theory given by (|l9|) and the classical one given by (18), since they both are in the confined 



phase. In particular, both theories should give a magnetic screening mass of order g'^T. 
It therefore seems likely that the classical theory is able to generate the approximatively 
correct non-perturbative behavior at the magnetic-mass scale, and that the replacement 
Oiw Oiw{T) in the classical expression (Q) yields a reasonable transition rate estimate in 
the full theory. 
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Appendix 



In this Appendix we briefly describe thermal properties of Chern-Simons number in the 
noncompact classical U(l) theory whose lattice Hamiltonian is 



(20) 



3,n 



where, in the usual lattice notation, j are sites and m, n are positive directions on a cubic 
lattice, and are lattice gauge potentials, whose conjugate momenta Ej^n are the electric 
fields. Note that upon substitution Uj^n = 6xp(zcj"0^„) and linearization ( [To| ) reduces to 
([20|) summed over color indices a. Imposing the Coulomb gauge 



-']~n,n) 



0, 



we eliminate longitudinal degrees of freedom (the electric field is transversal by the Gauss' 
law). In momentum representation (Bol) becomes 



H 



1 



E 



£p£-p + 4 2^ sm ( — 1 apa_ 



where the Cartesian components of the momentum are Pm. = '^'^Im/Li with integer 1 — 
Ll/2 < Im < Ll/2 on an L| lattice with periodic boundary conditions, a'_„ = (al)* is a 
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photon field with momentum p and polarization i, and £p = (<f!_p)* are electric fields in 
momentum space. It is easy to see that a mode with momentum p has an eigenfrequency 



In the same notation the Chern-Simons number is 

-^cs — 7^ n^lmn / , {(^j—n,n ~l~ ^j,n) ifi j~l+m,l ~\~ (^j+m,l Gj—m,l ^7— rn,0 

ozvr 

^eimn E e^(P'-P")/2 eos (^) cos (^) sm{pm)e\p)ne^ {-p)ia;alp, (21) 



p 



where the transverse polarization vectors e*(p) on the lattice have the properties 



/ \* sinf— ^sinf— ^ 

{-P) = (e\p)) ; e\pU^{-p)n = S'^; eHp)ie\-p)n = Sm - e^^-^") , ,l\ ^ 

V / ^,,sm (f) 



Using these properties, and performing a straightforward Gaussian integration with the 
thermal weight exp(—(3H), one obtains the autocorrelation function 



1 cos^ {uJpt) 

^P 

This expression can be explicitly evaluated for any finite lattice. The thermal width v 



of Abelian A'^cs is given by (22) at t = 0. In the limit L^, — > oo the sum in ( P2|) may be 
replaced by an integral, and becomes an extensive quantity: ~ 6.930 x 10~^/3~^L|. 
We used this estimate, together with the conversion rule (|l^), to obtain (0). Note that the 
long-time behavior of ( p2| ) depends on the order in which the thermodynamic Li oo 
and the t — > oo limits are taken: if the latter is taken first (which is the case of practical 
interest here) on a finite-size lattice, ( p2| ) will at long times fluctuate around u and not 
approach zero, as can be verified by taking time average of (^) over a long time interval. 

Another useful Abelian quantity is the thermal width of dfNcs- Taking time derivative 
of (21) and averaging over the thermal ensemble yields 

(*i>,=^i:n-^(|)- Jy. m 

where in the last step the large-volume limit was taken. 



References 

[1] G. 't Hooft, Phys.Rev.Lett. 37 (1976) 8; Phys.Rev D14 (1976) 3432. 

[2] F. Klinkhamer and N. Manton, Phys.Rev. D30 (1984) 2212. 

[3] D. Yu. Grigoriev and V.A. Rubakov, Nucl.Phys. B299, (1988) 671. 

D. Yu. Grigoriev, V.A. Rubakov and M.E. Shaposhnikov, Nucl.Phys. B326 (1989) 
737; Phys.Lett. B216 (1989) 172. 

[4] S.Yu. Khlebnikov and M.E. Shaposhnikov, Nucl.Phys. B308 (1988) 885. 



11 



p. Arnold and L. McLerran, Phys.Rev D36 (1987) 581. 

D. Bodeker, L. McLerran and A. Smilga, Really Computing Non-perturbative Real 
Time Correlation Functions, hep-th/9504123| . 



J. Ambj0rn, T. Askgaard, H. Porter and M. Shaposhnikov, Nucl.Phys.B353 (1991) 
346; Phys.Lett. B244 (1990) 479. 

P. de Forcrand, A. Krasnitz, and R. Potting, Phys.Rev. D50 (1994) R6054. 

A. Krasnitz, Thermalization algorithms for classical gauge theories, NBI report NBI- 



HE-95-25, |hep-lat/950702q , to be published in Nucl.Phys. B; Nucl.Phys. B (Proc. 



Suppl.) 42 (1995) 885. 

J. Ambj0rn and K. Farakos, Phys.Lett. B294 (1992) 248. 

J. Ambj0rn, K. Farakos, S. Hands, G. Koutsoumbas and G. Thorleifsson, Nucl.Phys. 
B425 (1994) 39. 

U.M. Heller, F. Karsch, and J. Rank, The gluon propagator at high temperature, 
Bielefeld report BLTP 95/21, |hep-lat/95050T^ . 

T. Reisz, Z.f.Phys. C53 (1992) 169. 

P. Lacock, D.E. Miller and T. Reisz, Nucl.Phys. B369 (1992) 501. 

L. Karkkainen, P. Lacock, B. Petersson, D.E. Miller and T. Reisz, Nucl.Phys. B395 

(1993) 733. 

K. Farakos, K. Kajantie, K. Rummukainen and M.E. Shaposhnikov, Nucl.Phys. B425 

(1994) 67; Nucl.Phys.B442 (1995) 317; Phys.Lett.B336 (1994) 494. 

F. Karsch and H.W. Wyld, Phys.Rev. D35 (1987) 2518. 

R. Berman, Thermal conduction in solids, Clarendon Press, Oxford (1976) 14. 



12 



